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Abstract 

Inelastic X-ray Scattering measurements of the dynamics structure factor 
have bene performed on vitreous Beryllium Fluoride (?;-BeF2) at T=297 K in 
the momentum transfer, Q, range (5=1.5-^10 nm^^. We find evidence of well 
defined high frequency acoustic modes. The energy position and linewidth of 
the excitations disperse with Q as oc Q and oc , respectively up to about one 
half of the first maximum of the static structure factor. Their magnitude com- 
pare favorably with low-frequency sound velocity and absorption data. The 
results indicate worth mentioning similarities of the high frequency collective 
dynamics of different network forming glasses such as U-B2O3 and u-Si02. 
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I. INTRODUCTION 



The interest in the provocative behavior of the glass-forming hquid BeF2 has been re- 
cently renewed in an effort to relate a viscosity anomaly to some water-like features, such 
as the existence of a negative thermal expansion coefficient region [|l|. Since the first ap- 
pearance of the Tg-scaled temperature dependence of the viscosity plot for glass-forming 
systems, it became clear that liquid BeF2 has an intriguing property: it shows a crossover 
between the two extreme strong and fragile behaviors. In particular, using existing viscosity 
data 1^ Angell and co-workers have shown that -in a T^-scaled plot- the high-temperature 
limit of BeF2 viscosity obtained by extrapolation should reach an un-physically small value 
Molecular dynamics simulations have recently been performed to explain the above 
observation in terms of a weak thermodynamic anomaly. In particular, a density maxi- 
mum has been predicted to occur at 2000 K and a density minimum at 1250 K, a 
temperature at which a considerable (~30%) rise in heat capacity was also predicted. 

Despite these intriguing thermodynamics anomalies, BeF2 is much less studied compared 
to other network forming glasses like the oxides Si02, Ge02, B2O3 etc., owing to the sub- 
tleties that it presents in its purification procedure (hygroscopicity, toxicity, corrosiveness 
etc.). 

From the experimental point of view the structure of BeF2 has been elucidated by x- 
ray ^ and neutron diffraction studies supporting the resemblance of BeF2 structure 
with that of Si02, i.e. a 3D network of corner sharing BeF4 tetrahedral units. The room 
temperature Raman spectrum of BeF2 has been measured by Galeener et al. where a 
comparison between this glass and other oxides is attempted. Finally, dynamic properties of 
BeF2 at low temperature and ultrasonic frequencies have been carried out |^ showing that 
i) between T=10 and 200 K, BeF2 has an acoustic absorption comparable to that of Si02, 
a) the temperature dependence of the sound speed of BeF2 (similar to that of Si02) shows 
anomalous behavior. In particular, in contrast to the monotonically decreasing sound speed 
with temperature rise which is usually found in glasses (B2O3, Ge02, Zn(P03)2), the sound 
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speed of BeF2 (and of Si02 as well) exhibits an initial drop followed by an upturn above 50 
K. 

In this paper, we present the first experimental determination of the acoustic properties of 
BeF2 in the Ti/z frequency range, by means of Inelastic X-ray Scattering (IXS). In particular, 
we have measured the dynamic structure factor, S{Q,uj), at room temperature (Tg=598 K) 
i) a.s a function of energy E for fixed values of Q in the range 1 -i- 8 nm^^; ii) as a function 
of Q for the fixed energies of and 7 meV. The obtained data allow us to extract some 
information about the high frequency sound dispersion and attenuation properties for the 
longitudinal acoustic excitations of this network forming glass. Specifically, we find that 
i) Well defined high frequency acoustic modes exist and propagate in the glass; ii) The 
excitation frequency Q{Q) disperses linearly with Q {fl{Q) = vQ); in) The excitation width 
T{Q) (FWHM) increases quadratically with Q, (T{Q) = DQ^); iv) The sound speed, 
and the sound energy absorption coefficient, a = 27cr/v, compare favorably with literature 
low- frequency ultrasonic data |^. Finally, a comparison of the sound absorption data of 
BeF2 with that of Si02 indicates striking similarities in the behaviour of the high frequency 
collective dynamics of these two systems. 

II. EXPERIMENTAL DETAILS 

Beryllium fluoride is a substance that presents considerable difficulties in its purification 
procedure. This stems from the fact that it is extremely hygroscopic, toxic and corrosive 
for conventional containers like fused silica tubes. For the above reasons the whole material 
handling operation took place in an inert atmosphere (nitrogen-fllled) glove box with a 
water content less than 2 ppm while the material was melted in gold-plated silica tubes of 
8 mm inner diameter. The BeF2 starting material was purchased from Alfa Aesar Chemical 
Co. with nominal purity 99.5 % which was not enough to obtaine a transparent glass/melt 
free of black spots. Thus, BeF2 was flrst sublimed under dynamic high vacuum and high 
temperature in graphite tubes. The material obtained after sublimation contained pure and 
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contaminated parts from the container. The proper amount of pure BeF2 was selected and 
placed into the gold-plated cylindrical silica crucibles with flat bottom. The crucibles with 
the pure BeF2 were placed in a quartz cell that was evacuated, partially filled with Ar, and 
flame sealed. This cell was then transferred into a furnace and heated up to the softening 
point of the glass (1100 K) for 1-2 hours to obtain complete homogenization. The melt 
was cooled down to room temperature in a controlled way in order to obtain a transparent 
glass, free of internal stresses. With this procedure we were able to obtain cylindrical BeF2 
glass samples with a length of about 1 cm, comparable with the absorption length for 21.7 
keV x-rays. Characteristic temperatures of BeF2 are: glass transition temperature Tg = 598 
K and melting point of both the quartz-like and cristobalite-like crystalline modifications 
Tm = 828 K. 

Data for the BeF2 glass have been collected at room temperature (T=297 K) at the 
IXS beamline ID28 of ESRF P,^. The experiment has been performed at fixed exchanged 
wavevector over a Q— region 1.5 ^ 10 nm^^ with a resolution (FWHM) of 6Q ~ 0.35 nm~^. 
The overall energy resolution (FWHM) has been set to 6E = 1.5 meV utilizing the (11 
11 11) reflection for the Si monochromator and crystal analyzers. A five- analyzers bench, 
operating in horizontal scattering geometry, allowed us to collect simultaneously spectra at 
five different values of Q. Each energy scan (—35 < E < 35 meV) took about 300 minutes, 
and several scans have 

been accumulated for a total integration time of about 500 seconds/point. Measurements 
at constant energy have also been conducted, scanning over the scattering angle. In this 
case the energy resolution was the same as in the fixed Q scans, while the Q resolution was 
increased, 6Q ~ 0.17 nm^^, by using tunable-width slit placed on the scattered beam path. 

III. DATA PRESENTATION AND DISCUSSION 

Fig. |1] illustrates a selection of experimental spectra accumulated at the indicated Q 
values as a function of the exchanged energy E. Since the incident flux on the sample slightly 
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varies during the acquisition of each scan, the data have been normahzed to a monitor signal 
for each frequency value, and then multiplied by the average monitor. The data of Fig. 1 
bear a close resemblance with those collected from other network forming glasses, such as 
Si02 [[l^iOl ^"^^ B2O3 |T^. A strong elastic peaks dominates the spectra, and the inelastic 



features appear only as weak shoulders on the tail of the resolution-broadened elastic line. 
In order to extract quantitative information on the excitations giving rise to the inelastic 
signal, we fitted the the experimental data with a model function for the S{Q, E) convoluted 
with the instrumental resolution. 

In fact, the actual experimental intensity, I{Q,E), is proportional to the convolution 
of the dynamic structure factor, Sq{Q,uj), {u = E/h), with the instrumental resolution 
function R{E). 

I{Q, Lu) = e{Q) J dcu'SgiQ, u')R{uj - u') (1) 

where e{Q) contains the efficiency of the analyzers, the atomic form factors and other 
angular-dependent correction factors. The true, quantum, Sq{Q,uj) can be approximately 
related to its classic counterpart by 

SQiQ^^)^Y^^^^S{QM (2) 

An useful expression for the (classical) dynamic structure factor can be obtained by 
recalling that its Fourier transform in the frequency domain, i.e. the density fluctuation 
correlation function E{Q,t) = {6pQ(t)6pQ{0)), obeys a generalized Langevin equation |13[: 

F{Q,t)+ujlF{Q,t)+ I m{Q,t-t')F{Q,t')dt' = (3) 

Jo 



where Uo is a parameter related to the static structure factor S{Q) ||13[, and ■m{Q,t) is the 
"memory function". In this "exact" expression all the difficulties associated to the calcula- 
tion of F{Q,t) have been transferred to the determination of m{Q,t), with the advantage 
that the first two sum rules for S{Q,uj) are automatically satisfied [|T3|. By Fourier trans- 
forming Eq. (H), it is easy to show that: 



S{Q,uj) 1 oo^m'iQ^oo) 



(4) 



where m{Q,u!)=m'{Q,u!)+i'm"{Q,uj) is the time Fourier transform of the memory function 
m{Q,t). In the cjTq, >>1 hmit, a hmit that is certainly vahd in the present case of a glassy 
sample, the memory function can be approximated by the sum of a constant term, A'^[Q), 
-which reflects the frozen a-process- plus a function showing a very fast decay at short 
times. The latter contribution to the memory function -often referred to as "microscopic" 
or "instantaneous"- is usually represented as a delta-function with area 2r{Q). Therefore: 

m{Q,t) = 2T{Q)6{t) + AUQ) (5) 

and hence Eq. reads as: 

'f ^ ^ i-/q ^\QnQ) 1 



S{Q) 



where n(Q) = A/Af(QyT~cuf, and fQ=l-uj1/VL^{Q) is the non-ergodicity factor. The expres- 
sion m Eq. I is the sum of an elastic line (the frozen a process) accounting for a fraction 
fq of the total intensity, and of an inelastic feature which is formally identical to a Damped 
Harmonic Oscillator (DHO) function. The parameter Vt{Q) coincides with the maximum of 
the longitudinal current correlation function, J{Q,uj) = {u"^ / Q'^) S {Q , u) , and it is related 
to the apparent sound speed value q(Q) = Q{Q)/Q, while T{Q), the excitation width, is 
related through its low Q value to the acoustic absorption coefficient a = 2* 7rr(0)/Q(0). 
The model function, I^\Q,uj) utilized to fit the experimental data will be, therefore 

/*^(g, u) = aiQ) J j^^^Riuj - X (7) 



fQ6{uj') + ^-ir^(^/2_f^^)2+^/2r^ 



du' 



which, apart from an intensity factor a{Q) = e{Q)S{Q), contains three free shape parame- 
ters: Qq, Tq and fq. 

In Fig. |l|, together with the raw data (circles), we report the the elastic contribution 
(dashed lines), the fits to the data (full lines) and the inelastic contribution (dotted lines). 
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This latter contribution clearly shows a dispersing behavior despite the considerable sound 
attenuation {T{Q*) = n{Q*) at Q* ^ 3 nm'^). 

In Fig. 1^ (upper panel) we report the best fit values of Q{Q) (dots) along with the 
dispersion obtained in the low frequency (ultrasound) limit (dotted line, corresponding to 
the adiabatic sound speed Cq = 4570 m/s |T^). To represent the Q-dependence of fl{Q) 
we utilized the expression Q{Q) = ci{0)Qosin{Q / Qo) in order to keep into account for the 
bending of the dispersion curve associated to the existence of a quasi-Brillouin zone. This 
procedure yielded for the sound speed q = 4790 ±50, a value that, as in the case of vitreous 
silica |Tl[, exceeds its low frequency value. The higher value of the high frequency sound 
velocity could be the signature of the crossing of some relaxation process. Such a process 
can not be the structural relaxation, since in the glass it lies in the snh-Hz region, while it 
can be ascribed to the action of a microscopic process -produced by the topological disorder 
that is present in the glassy phase, as demonstrated in the case of simulated monatomic 



glasses []T5|,|T6[. This hypothesis implies that the assumption of an instantaneous decay of 
the memory function in eq. |^ may be inadequate; unfortunately, however, the low magnitude 
of the inelastic to elastic ratio does not allow us for a more realistic ansatz for the memory 
function shape. The excitation width Tq is reported as a function of Q in the lower panel of 
Fig. H. Its Q-dependence is compatible with the law already found in all the glasses where 
the dynamic structure factor has been measured in the mesoscopic region: Fg = DQ^. 
Actually, a fit to the data with the law Tq = DQ'^ (full line) results in 7=1.96±0.05 and 
D=1.05±0.05 meV/nm-2. 

In order to further illuminate the previous results, the dynamic structure factor has also 
been measured at a fixed exchanged energy value as a function of the momentum transfer. In 
the constant Q scans (as those reported in Fig. 1), indeed, the inelastic signal always appears 
as a shoulder of the resolution broadened central line, so that the choice of an appropriate 
model for the S{Q,uj) can be, in principle, neither easy nor unique [|T^. In a constant E 
scan, on the contrary, the inelastic signal is more clear, since the (resolution broadened) 
elastic contribution appears as an almost fiat, Q- independent, background. Moreover, the 



parameters determined trough the best fit to the energy scan data can be utihzed to build- 
up a curve that can be compared to the constant E experimental data |T^. In Fig. ^ we 



present, as an example, the result of such a comparison for a scan at E=7 meV: in the upper 
panel we show the raw experimentally measured elastic, I{Q,E = meV), and inelastic, 
I{Q,E = 7 meV) data. Obviously, the latter data contain an elastic contribution as a 
consequence of the finite energy resolution. In order to subtract this elastic contribution we 
used a series of spectra at constant Q taken between Q=23 and 32 nm~^. After aligning and 
scaling the experimentally determined resolution function to the elastic peaks in the spectra 
at high Q, we estimate the relative intensity between the elastic and inelastic signals at the 
energy transfers utilized in the constant-i? spectra. The elastic to inelastic intensity ratio 
has been obtained at Q=13, 16, 19, and 22 nm~^. These ratios allow us to put in the correct 
relative scale the spectra taken at E=7 meV and at E=0. This normalization procedure 
is used to derive the inelastic part of the S{Q, E) by the subtraction of the normalized 
elastic contribution from the total scattered intensity. The difference spectrum at E7 meV 
is reported in the lower panel of Fig. ^ (circles) together with the error bars as derived from 
the counting statistics. In this spectrum the existence of a defined Brillouin peak is highly 
emphasized. The dashed line in the lower panel of Fig. [Represents the function predicted on 
the basis of Eq. |^, with VLq = ci{0)Qosin{Q / Qo) and Tq = DQ'^ (as derived from Fig. H), 
where the only adjustable parameter is now an intensity factor since q(0), Qq, D and 7 are 
determined from the fit results of the energy scans. Although the peak position turns out to 
be slightly underestimated, the proposed fitting model appears to capture the main features 
of the experimental data. 

The convergent results obtained by analyzing the two independent data sets, i. e. the 
constant-Q and -E scans, allow us to establish the appropriateness of the approximations 
introduced in the memory function of the density fluctuations, and, more specifically, of the 
results reported in Fig. ^ for VLq and Fg. 

One of the most intriguing topics, as far as the glassy dynamics is concerned, is the 
frequency (or Q) dependence of the sound attenuation. Despite that in the IXS window 
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(i. e. at Q values ranging from ^1 to ~10 nm~^, or energies ranging from ^5 to ~50 meV) 
a law describes the excitation broadening, at lower frequencies the situation becomes less 
clear due to the presence of relaxation processes. An interesting comparison between BeF2 
and Si02 is presented in Fig. |^ where we report the T{Q) parameter for these two glasses 
over a wide range of frequencies as derived from IXS and other literature data. Both glasses 
exhibit a behavior above 35 GHz, while below this value a Q^-^ power law well represents 
the available experimental data for Si02 . 

Finally, in Fig. H, we present the non ergodicity factor as determined by IXS in vitrous 
sihca and Beryllium Fluoride. Indicating with lei = S{Q)f{Q) and linei = S{Q)[1 — f{Q)] 
the integrated elastic and inelastic spectral contribution (determined by the fit), the non 
ergodicity factor reads as f(Q) = T^^irr- As shown in the figure, the f{Q) of BeF2 is 
in good agreement with the Si02 data when reported as a function of T/Tg. The results 
reported in Figs. ^ and § imlpy the existence of an interesting analogy between these two 
network forming glasses as far as their high frequency dynamics is concerned. 



IV. CONCLUSIONS 

In conclusion, a room temperature IXS study of glassy BeF2 has been undertaken in the 
present work. In accordance with all previous studies in strong and fragile glasses, evidence 
has been presented for well-defined propagating (high frequency) acoustic modes, whose 
frequency position and linewidth scale as oc Q and oc Q^, respectively. The longitudinal 
speed of sound for glassy BeF2 has been estimated to exceed its low frequency (ultrasonic) 
limit by almost 5%; a case analogous to that found in studies of vitreous silica. The ex- 
trapolation of the high-frequency linewidth conforms nicely with the value obtained from 
ultrasonic studies, and exhibits a scenario similar to that of vitreous silica. Another similar- 
ity which deserves further study is that the temperature dependence (in a T^-scaled plot) of 
the non-ergodicity factor, as determined from the ratio of the elastic to the total scattered 
intensity, follows the behavior exhibited by Si02 (Fig.^ while for less strong glasses the 
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drop of f{Q) is much faster with increasing temperature. Unfortunately, for BeF2 it is up 
to now available only one point in such a T^-scaled plot, which however coincides with the 
f{Q) data for silica. Further temperature-dependence studies on BeF2 are expected to shed 
more light on this issue. 
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FIGURES 
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FIG. 1. Selected IXS energy spectra of the BeF2 glass at room temperature and fixed values 
of the exchanged wavevector (open dots o). The solid lines represent the best fits with Eq. 0, 
while dashed and dotted lines account for the elastic and inelastic contributions, respectively, as 
determined by the fitting procedure. 



Fig. 1 - T. Scopigno et al, "High frequency acoustic modes..." 
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Q(nm') 

FIG. 2. Q— dependence of the sound frequency and attenuation. Upper panel: experimental 
sound dispersion (full dots) and best fit through a sine function (full line). The adiabatic sound 
speed, measured by ultrasounds, is also reported (dotted line). Lower panel: Acoustic attenuation 
as determined by the Brillouin linewidth (full dots). The solid line is a power law fit to the 
experimental data and gives Tdho = 1.05Q^'^^ 



Fig. 2 - T. Scopigno et al, "High frequency acoustic modes... 
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FIG. 3. Fixed energy IXS Q— scan. Upper panel: inelastic {I{Q,E = 7 meV), full dots) 
and elastic {I{Q, E = meV), open dots). Lower panel: genuine inelastic signal determined as the 
difference of the above reported scans (normalized according to Ref. fl^), compared to a theoretical 
DHO with parameters determined as discussed in the text. 



Fig. 3 - T. Scopigno et al, "High frequency acoustic modes..." 
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FIG. 4. Sound attenuation of BeF2 and Si02 determined with different techniques (the right 
scale is for Si02 data only). The full lines are and Q^'^ power laws. The inset details the IXS 
frequency window 



Fig. 4 - T. Scopigno et al, "High frequency acoustic modes... 
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FIG. 5. Non-ergodicity factor of BeF2 and Si02 in a Tg scaled plot. The reported data have 
been determined for momentum transfers ranging from 1.6 to 2.5 nm~^, a region where f{Q) turns 
out to be rather constant. The full line is a guideline to the eyes only. 



Fig. 5 - T. Scopigno et al, "High frequency acoustic modes... 
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